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The asymmetry in the forward region production cross section of 
is calculated using the heavy quark recombination mechanism for pp col¬ 
lisions at 7 TeV. By suitable choices of four nonperturbative parame¬ 
ters, our calculated results can reproduce those obtained at LHCb. We 
hud Ap ~ —1% when integrated over 2.0 GeV < pr < 18 GeV and 
2.2 < T] < 4.75, which agrees with Ap = —0.96 ±0.26 ±0.18% as measured 
by LHGb. Furthermore, the calculated distributions in rj and pt agree 
reasonably well with those obtained at LHGb. Using the heavy quark 
recombination mechanism, we also make predictions on the production 
asymmetries of /A“ and A^A^ for pp collisions at 7 TeV and 14 TeV 
in the forward region. We hnd that the integrated asymmetries for these 
A baryons in the LHGb region are of the order of ~ 1 — 3% and should 
be measurable. 
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1 Introduction 


One of the simplest signals for CP violation in charm is obtained by comparing 
partial decay widths of charm mesons to those of anticharm mesons. While CPT 
symmetry requires the total widths of D and D to be the same, the partial decay 
widths P(D —?• /) and P(i5 —)■ /) are different in the presence of CP violation, which 
is signaled by a nonzero value of the asymmetry 


riD^n-TiD^f) 

r{D ^f) + r(D^f) 

This signal is reasonably robust for /D~ mesons, provided that the number of de¬ 
caying particles and antiparticles is the same. However, at the Large Hadron Collider 
(LHC), the number of produced and D~ mesons might not be the same due to 
the fact that the initial state contains two protons. With CP-violating asymmetries 
expected to be at the per mille levels pQ, it is important to examine the production 
asymmetry of D mesons both experimentally and theoretically. 

Indeed, fixed-target experiments have already observed large asymmetries of charmed 
mesons and baryons in the forward region. In hadroproduction, the charmed hadrons 
are preferentially produced with a light valence quark of the same type as what ap¬ 
pears in the hadronic beam, for example |2]. This has been termed the “leading 
particle effect”. More recently, a similar asymmetry in production, defined as 


_ a{D+)-a{D-) 
a{D+) + a{D~) ’ 


( 2 ) 


has been measured in the forward region to be —1% by the LHCb Collaboration [3]. 
In Section we report our hndings of explaining this production asymmetry at 
LHCb by the heavy quark recombination mechanism |1]. In Section using the 
heavy quark recombination mechanism, we predict the production asymmetries of 
A+/A“ and A^/A^ for pp collisions at 7 TeV and 14 TeV in the forward region [S]. 
We conclude in Section IH 


2 production asymmetry from heavy quark re¬ 

combination mechanism 

Factorization theorems of perturbative QCD [6] state that heavy hadron production 
cross section can be written in a factorized form. At the LHC, the cross section for 
producing a D (cq) meson in a pp collision, at leading order in a l/pr expansion, is 
given by 

da[pp -P- D + X] =fi/p ® fj/p (g) da[ij -)■ c + X] g Dc-^d , (3) 

ij 
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where fi/p is the parton distribution function for parton i in the proton, da{ij —)■ 
c+X) is the partonic cross section and -Dc-s>d is the fragmentation function describing 
hadronization of a c quark into a D meson. The corresponding equation for D is 
obtained by replacing c by c and D hj D. Owing to charge conjugation symmetry 
and that fc/p = fc/p, perturbative QCD Eq. ([^ predicts that Ap = 0, which is at 
least true at leading order in the 1/pr expansion. 

To reconcile the experimental observations with QCD, we note that there are 
corrections to Eq. (|^ that scale as powers of Aqcd/’^c and Aqcd/pt- In princi¬ 
ple, one can expect nonvanishing power-suppressed contributions to Ap at low pt- 
A QCD-based model for these power corrections is the heavy quark recombination 
mechanism miHiEiE]. In this scenario, a light quark involved in the hard scattering 
process combines with the heavy quark produced in that interaction to form the final 
state meson, leading to corrections of order Aqcd^c/Pt- This contribution to the 
cross section is given by (Fig. UKa)) 

da[D] = dd[qg -> (cq)" + c]/3|(c())” -> £>], (4) 

where (eg)"’ indicates that the light quark of flavor q with momentum of order Aqcd 
in the c rest frame is produced in the state n, where n labels the color and angular 
momentum quantum numbers of the quark pair. The cross section is factored into 
a perturbatively calculable piece dd[qg —)■ (eg)"' -|- c] and a nonperturbative factor 
p[(cg)" —?• D] encoding the probability for the quark pair with quantum number n 
to hadronize into a hnal state including a D. The perturbative piece was calculated 
to lowest order in [7j. Equation (|^ must then be convoluted with the proton parton 
distribution functions to get the hnal hadronic cross section. Besides the qg —)■ 
(eg)" -|- c process, there are also contributions from gc —)■ (eg)" -|- g, as shown in 
Fig. (b). Using the method introduced in [7], we calculate the partonic cross 
sections from initial state charm jl]. 

The c quark in Eq. (|^ could fragment into a D meson. Thus, to get the full 
rate due to recombination for producing D mesons, we also need to account for the 
contribution where a light antiquark comes from the proton, while the c fragments 
into a D. We thus have three contributions: 

(а) da[D] = dd[qg —)■ (eg)" -|- c]p[(cg)" —)■ D ], (5a) 

(б) dd[D] = dd[qc (eg)" + g]p[(cg)" ^ D ], (5b) 

(c) da[D] = dd[qg (eg)" + c]p[(cg)" ^ H] ® , (5c) 

where H can be any hadron. The recombination cross section for producing a D 
is obtained by taking the charge conjugate of the above equations. Below, we will 
neglect C violation and take p[(cg)" —)■ D] = p[(cg)" —)■ D]. For simplicity, in process 
(c) we will restrict hf to be D or D* only and sum over q = u,d and s with SU{3) 
flavor symmetry assumed. 
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Figure 1: Diagrams for production of a D meson by the heavy quark recombination 
mechanism for (a) qg —)■ (eg)” + c and (b) gc —)■ (eg)” + g. Each process has hve 
diagrams. Single lines represent light quarks, double lines heavy quarks, and the 
shaded blob the D meson. 


The leading contributions to productions of mesons by heavy quark recombi¬ 
nation consists of four possible options of n\ 

p‘r = p[cd(‘s“) ^ D*] , p\' = p|ca(=sW) ^ £>+] . 

rf” = />[cS(‘S®) ^ D *], Pi' = plcd^sf) ^ £>+]. (6) 

These nonperturbative parameters must be extracted from data. Neglecting and 
, the combination pf^+pf^/8 was determined to be 0.15 by htting to the E687 and 
E691 hxed-target photoproduction data [8]. Neglecting pg™', and pl^, the parame¬ 
ter pf^ was determined to be 0.06 by htting to data from the E791 experiment |9]. In 
this paper, we take p^”* ~ 0.06 and pg™ ~ 0.7. It turns out that these two contribu¬ 
tions only account for ~ 30% of the measured asymmetry Ap = (—0.96±0.26±0.18)% 
at LHCb in Ref. [3]. Therefore, we include pl^ and pg'^ and choose values of similar 
size as the spin-matched parameters. We also include feed down from D*^. From 
heavy quark spin symmetry, we have 

p[cd(^4"^) ^ D+] = p[cdCs{^^) D *+], 

p[cd(3Ri"^) ^ D+] = p[cd(^Rf^) ^ D*+]. (7) 


We use MSTW 2008 LO central PDFs with rric = 1.275 GeV and the Peterson 
parametrization for the fragmentation function is used for Dc-^h'- 


Dc-^h{.z) 



Nh 


( 8 ) 


Cc ~ {'mq/mc)‘^ was measured to be 0.062 ± 0.007 for the D*^ meson m- Charge 
conjugation symmetry and approximate heavy quark symmetry implies that Cc is 
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(a) 



(b) 


Figure 2: Asymmetry in production Ap as a function of (a) pseudorapidity r] 
and (b) transverse momentum px in 7 TeV pp collisions. The data points are from 
LHCb [3]. The grey band is obtained by varying the p’s in the intervals 0.055 < 
pj”^ < 0.065, 0.65 < pI'" < 0.8, 0.24 < p^'^ < 0.30 and 0.24 < pg-^ < 0.30 respectively. 
The dashed lines are from varying 0.055 < Cc < 0.069. 


approximately the same for and D*^. We will take Cc = 0.06 for both and 
D*^. Nh are determined by the averages of the measured fragmentation probabilities 
listed in [T2]. For the perturbative QCD rate, Eq. ([^, which has no asymmetry if we 
ignore C violation but enters into the denominator of Eq. ([^, we use the LO cross 
section. The factorization scale is set to be p/ = + m^. 

When integrated over 2 GeV < pr < IS GeV and 2.2 < p < 4.75, excluding the 
region with 2 GeV < pt < 3.2 GeV, 2.2 < p < 2.8, the asymmetry Ap for is 
found to be -0.88% < Ap < -1.07% with 0.055 < < 0.065, 0.65 < p|™ < 0.8, 

0.24 < < 0.30 and 0.24 < pg-^ < 0.30. Figure 0 shows Ap as a function of 

pseudorapidity p and transverse momentum p^ of the 7W mesons as predicted by the 
heavy quark recombination mechanism. Data from Ref. |3] are shown as well. The 
grey band is from varying the p parameters within the ranges above. The dashed line 
is obtained using the central value of the p parameters and varying ec within its error 
bars. The calculated distributions are reasonably consistent with the data. Figure 
shows the independence of Ap on the factorization scale pf when pj is varied the the 
range ^\/pt + < l^f < ‘^Vpt + The scale dependence is significant (~ 100%) 

at the high p and low pt ends. An NLO calculation for the heavy quark recombination 
mechanism in the future will presumably reduce this theoretical uncertainty. 
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Figure 3: Asymmetry in production Ap as a function of (a) pseudorapidity r] 
and (b) transverse momentum pt in 7 TeV pp collisions. The data points are from 
LHCb [3]. The grey band is obtained by varying /ij in the interval \\/pt + '^1 < 
Pf < '^\/pt + The p’s are taken to be central values of the ranges in Fig. i Cc is 
taken to be 0.06. 


3 A+/A^ and A2 /a|J production asymmetry from heavy 

quark recombination mechanism 

Similarly, one can apply the heavy quark recombination mechanism to predict pro¬ 
duction asymmetries for baryons. The contributions to production of Aq (udQ) from 


recombination are given by 


(a) 

ddlAq] = da[qg (Qq)^ + Q]p[(Qg)’" ^ Aq] , 

(9a) 

(b) 

daiAq] = da[Qq (Qg)" + g]v[{QqT ^ Aq] , 

(9b) 

(c) 

da[Aq] = ^ da[qg (Qg)*" + Q] Hmeson] 



^ ^ meson 

(9c) 

(d) 

da[Aq\ = ^ da\pg {QqY + Q] ^ vHQqT Hbaryon 

1 ® ■ 


^ baryon 

(9d) 


Here we include both mesons and baryons in the opposite-side recombination (c) and 
(d). For simplicity, we will take 77 to be a low-lying heavy hadron. Thus, for A^ 
production we will take Hmeson to be either 77 or 77*, and Hf,aryon be any baryon from 
the lowest mass = \ and | heavy baryon SU{3) flavor multiplets, and similarly 
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for Afe production. We will also assume SU{3) flavor symmetry. All the p’s and p’s 
scale as Aqqd/^q- The leading nonperturbative parameters ri[{Qq)"‘ —)■ Aq] are 

h3 = v[Q(lCS^^^) Aq] , fjs = v[QqCs[^^) Aq] , 

Ve = Aq] , fje = v[QqCs[^^) Aq] . (10) 

Contributions of feeddown from heavier baryons in processes (a) and (b) can be taken 
into account by using the inclusive parameter rjinc- 

Vinc[iQqT —t Aq] = p[(Qg)” —)■ Aq]+ ^2 VHQqT Ht}aryorL\B^Hfyf2^yQ^ ^ Aq-|-^] 

^baryon^^Q 

( 11 ) 

Here again we assume Hbaryon is a member of the lowest mass = 2 i heavy 

baryon SU{3) flavor multiplets. 

For the four p’s for D, we use the ranges in Section Best single-parameter £t 
to A^ asymmetry in hxed target experiments gives = 0.058 for A^ [TU]. We 

will take rjs^inc = r]e,inc = he,me = 0 and 0.052 < fjs^inc < 0.064 for Ac. For the ps for 
Ab and ps for B, we simply multiply the Ac and D counterparts by the theoretical 
scaling factor rric/mb. We use MSTW 2008 LO central PDFs with rric = 1.275 GeV 
and rrib = 4.18 GeV. The fragmentation function Dq^p^^ is taken as 

=/aq5(1-^), (12) 

where /aq is the inclusive fragmentation probability. This form of fragmentation 
function was found to be better than the Peterson form when fitting to hxed target 
A+ /A~ asymmetry data [lO]. We take fp+ = 0.101, which is the average of the values 
listed in [12]. f>.o is taken to be 0.09 from |T3j. The factorization scale is set to 

be ,, = Figueeege.d§eh™ the tapidtty attd ttanevetee —un. 

distributions of Ap for A+/A^ and A° / A^ in pp collisions at 7 TeV and 14 TeV 
respectively in the forward region. The asymmetry is signiheant at the high-rapidity 
and low-pT ends (~ 2 — 15%). 

Shown in Fig. are the GMS data na, which are the rapidity and transverse 

_Q 

momentum distributions of ct(A^)/(t(A|]) for 7 TeV pp collisions in the kinematic 
region 0 < p < 2 and 10 GeV < pr < 50 GeV. The GMS data, despite of the 
large error bars, do have a slight trend of surplus of A° over A^ in the regions of 
high rapidity and low transverse momentum respectively. However, with the values 
of p’s we used above, the asymmetry predicted from the heavy quark recombination 
mechanism is negligible in this kinematic region. To see the limit of the heavy quark 
recombination mechanism in explaining the data, in Fig. |^we also plot the prediction 
with larger values of p’s. Here all pme’s are set equal to each other, with the range 
being AQcnj'^h ~ 0.2 < pj„c < 1- The ranges of the ps are as those used in Fig. 
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Although the prediction shows a significant asymmetry ~ 10% at the high-rapidity 
and low-pr ends, it still fails to hit the bin with the largest rapidity. We believe that 
our previous predictions with smaller rjincS (Figs.|4||^ are more reasonable since those 
values of r^’s were obtained from better fit to fixed-target experiments. Moreover, the 
surplus of Af, over A° at ?/ ~ 1 and pr > 20 GeV in the CMS data cannot be explained 
by any existing model. We hope that data from LHCb in the future will settle the 
issue. 


4 Conclusion 

In summary, we have calculated the asymmetry using the heavy quark recombina¬ 
tion mechanism for production at the LHCb experiment. The measured asymmetry 
of Ap = —0.96 ± 0.26 ± 0.18% in the kinematic range 2.0 GeV < pr < 18 GeV and 
2.2 < 1 ] < 4.75 |3j can be reproduced using reasonably sized nonperturbative param¬ 
eters Further, the px and rj distributions are simultaneously reproduced by 

the heavy quark recombination mechanism. 

We have also used the heavy quark recombination mechanism to calculate the 
production asymmetries for A+ /A~ and A^/A^, at the LHCb experiment. The differ¬ 
ential distributions are significant at the high-rapidity and low-pT ends (~ 2 — 15%). 
The integrated asymmetries in the LHCb region are of the order of ~ 1 — 3% and 
should be measurable. 
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